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ABSTRACT 

The classification of galaxy clusters according to their X-ray appearance is a powerful tool to discriminate 
between regular clusters (associated to relaxed objects) and disturbed ones (linked to dynamically active systems). 
The compilation of the two subsamples is a necessary step both for cosmological studies - oriented towards 
spherical and virialized systems- and for astrophysical investigations - focused on phenomena typically present in 
highly disturbed galaxy clusters such as turbulence, particle re-acceleration, magneto-astrophysics . 

In this paper, we review several morphological parameters: asymmetry and fluctuation of the X-ray surface 
brightness, hardness ratios, X-ray surface-brightness concentration, centroid shift, and third-order power ratio. 
We test them against 60 Chandra-like images obtained from hydrodynamical simulations through the X-ray Map 
Simulator and visually classified as regular and disturbed. 

The best performances are registered when the parameters are computed using the largest possible region (either 
within /?5oo or 1000 kpc). The best indicators are the third-order power ratio, the asymmetry parameter, and the 
X-ray-surface-brightness concentration. All their combinations offer an efficient way to distinguish between the 
two morphological classes achieving values of purity extremely close to 1 . A new parameter, M, is defined. It 
combines the strengths of the aforementioned indicators and, therefore, resulted to be the most effective parameter 
analyzed. 

Subject headings: clusters: general — clusters: ICM — clusters: cosmology — X-rays: clusters 



1. INTRODUCTION 

For over two decades, the X-ray-galaxy-cluster community 
has been searching for methods to suitably quantify the dynam- 
ical status of clusters based on their X-ray morphology ( |Jones| 
|& Forman|1992| l. Initially, the morphological classification was 
mostly conducted by eye and strictly connected to the presence 
and characterization of substructures. For example, the mor- 
phological classes identified by Jones & Forman ( 1992| l were: 
'single', 'double', 'primary with small secondary', 'complex', 
'elliptical' (referred to the X-ray contours), 'off-center' (either 
presenting a difference between the centers in optical and X- 
ray or showing an X-ray tail extended only in one sector off the 
X-ray peak), and 'galaxy' (when the X-ray emission is dom- 
inated by the central galaxy emission). A very similar clas- 



sification was proposed by Baier et al. (T996JI who compared 
the substructures detected in X-ray with those known in op- 
tical wavelengths (see also Kolokotronis et al. 2001 1. These 



visual estimates appeared to be robust ( Jones & Forman|1 999) 



when compared with a more quantitative measurement of X 



ray-morphology disturbance such as the power ratios ( Buote & 
|Tsai|1995|[T996l l, or axial ratios ( |Mohr et al. 11993] ) . 

Since its beginning the detection of X-ray substructures has 
been associated with dynamically active systems (se e also[5 lezak 
eTaT1[T994l |Gomez et al.|[T997i |Rizza etaLl[T998l |Kolokotro-| 



nis et aL 2001 ; S chuecker et al.||2001| l. Indeed, as demon- 



strated through hydrodynamical simulations, a significant long 
time, spanning from 2Gyr to 4 Gyr, is required to re-establish 
a situation of equilibrium after a major merger (Lacey & Cole l 



1993 



et al. 



Roettiger et aL][T996l [19971 l Poole et al - 



2006; Nelson 



2012| >. However, even if the presence of substructures 



is certainly indicative of a dynamically active system (Rich- 
stone et al.|1992 1 the reverse is not always true. A cluster pre- 



senting either a strongly elliptical shape (Buote & Tsai 1996; 
Pinkney et al. 1996 |Phonis||2002l > or a variation of its X-ray 
centroid ( Mohr et al.|1995) is very likely an active cluster even 
if it has not recently interacted with a merging massive system 
and, therefore, does not present any substructure. Moreover, 



the visual classification becomes challenging in prospect of the 
large number of data provided by future survey. A need for 
more robust indicators that can quantify objectively even small 
deviations from a perfectly regular and spherically-symmetric 
emission is, therefore, highly desirable. 

Among the various methods proposed, we recall some based 
on wavelets analysis ( |Pierre & Starck|1998{|Slezak et al.|19 94), 



on the Minkowski junctionals (|Beisbart et al.|2001 1, on axial ra 



tios ( |Mohret al.|1993) , o n centroid-shift (e.g.JMohr et al.|1 995), 
and on power ratios (e.g. Buote & Tsai T995J. Several morpho- 
logical indicators, at the beginning, were used to test cosmolog- 
ical model through a comparison of the X-ray morphological 
measurements performed in X-ray cluster surveys and in hy- 
drodynamical simulated samples (Evrard et al. 1993 Jing et al. 



1995;,|Dutta|1995HMohr et al.|1995HTsai & Buote|1996HCrone 



et al.[l996HBuote & Xu|1997HThomas et al.|1998MBuote|1998t 
Valdarnini et al. 1999; Suwa et al. 2003; Ho et al. 2006). Their 



usage in this context, has been significantly reduced recently 
as they are not as robust as expected for the level of precision 
required by current and up-coming cosmologically-driven sur- 
veys. However, a way to classify the X-ray appearance is still 
of great interest to generate subsamples of either very disturbed 
or very relaxed objects. 

The first category provides the fundamental material for as- 
trophysical investigations of many microphysical processes. Merg 



ers are activating turbulence in the intracluster medium ( |Hall 



man & Jeltema 2011 Vazza et al. 2012 1, in particular in the 



sloshing core ( ZuHone et al.|2012a Roediger et al. 12012} . The 
detection and statistics of cold fronts formed from the sloshing 
gas are a direct test of presence of thermal conduction ( |ZuHone| 
et al. 20 1 2b) > . Studies of cold-front edges will allow also to 



estimate the level of viscosity of the fluid through detection 
(or no-detection) of Kelvin-Helmotz instabilities (Marke vitch| 
|& Vikhlinin|2007| >. Furthermore, merging clusters show corre- 
lation between the X-ray emission and the presence of giant ra- 
diohalos ( fBuote|2001||Cassano et al.|2010[|HaIlman & Jeltema| 
201 l| l. In this respect, the creation of highly-disturbed clus- 



ter sample will enable the examination of particle acceleration 
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mechanisms resulting in non-thermal radio emission (Mann & 
|Ebeling|2012l >. 

Regular clusters, instead, can be used to calibrate scaling re- 
lations. Since their X-ray mass derived by assuming hydro- 
static equilibrium is believed to be more robust ( Rasia et al.| 
2004] |2006l |Nagai et al.||2007j |Piffaretti & Valdarnini||2008| 
Lau et al.|2009 1 and closer to the gravitational lensing masses 



|Ma 



Meneghetti et . al.pOTOl |Rasia et al.||2012l |Chon et al.pOT? 



duced scatter in scaling relations (e.g. Rowley et al. 


2004| |Cros-| 


ton et al.||2008 Yang et al.||2009 Rasia et al. 


2011 


). Irregular 



systems, on the contrary, can be used in scaling relation only 
if their level of substructure is known (see recent method pro- 
posed by Andrade- Santos et al.|20T2 i. Indeed, even if typically 



irregular systems at fixed X-ray temperature tend to be less lu- 
minous ([Torri et al.||2004"t |Kay et al.||20"07l [Hashimoto et al.| 
|2007i [Ch 



on et al.|2012[ |Andrade-Santos et al.||2012|) and less 



massive ( O'Hara et al.|2006)|Valdarnin i 2006 ; Ventimiglia et al. 
|2008l[2012f|Chon et al.|2012||Andrade-Santos et al.|2012) , the 
parametrization of the degree of substructures can be used to 
correct their positions in the L—T or M— T planes (Ventimiglia 
|et al.|2008) . 

Morphological parameters will be therefore useful to catego- 
rize the large amount of data that will be collected from future 
surveys and to parametrize their X-ray morphology. With this 
perspective in mind, for the first time, we test several morpho- 
logical indicators against the same ensemble of synthetic X-ray 
data. The final goal is to identify the parameter(s) that most 
efficiently allow for the classification of regular and irregular 
objects. We establish these classes by visual analysis of the im- 
ages, i.e. following a pure observational approach. At first, we 
distinguish three classes: regular, disturbed, and 'mix'. In the 
last Section we subdivided the last category into semi-regular 
and semi-disturbed systems. 

The simulated sample and the X-ray synthetic catalogue are 
presented in Section 2, while the morphological parameters are 
introduced in Section 3. Each parameter is tested against the 
visual classes of regular and disturbed objects (Section 4). To 
allow for statistical evaluations, we extend the analysis to the 
whole sample when studying the improvement generated by 
combining two parameters (Section 5). Conclusions follow in 
Sections 6. In the appendix, we present the soft band images of 
all our clusters. 

2. SIMULATIONS 

The 60 X-ray images analyzed in this paper reproduce Chan- 
dra-like event files of 20 simulated clusters observed along three 
lines-of-sight each. This set of synthetic observations have been 
also used in Rasia et al. (2012, hereafter R12) and Rasia et al. 
(in prep.). Here we briefly describe the main characteristics of 
the simulations and of the X-ray catalogue. Further details can 
be found in the aforementioned publications. 

2.1. Simulated clusters 

The simulated clusters used in this work are obtained as fol- 
lows. Twentynine Lagrangian regions have been selected from 
a large (1 Gpc 3 /T 3 ) DM-only cosmological box. Twenty five of 
these regions are centered on the most massive objects of the 
parent simulation while the remaining are around group-size 
halos. The regions have been subsequently re-simulated with 
an increased resolution in mass following the Zoomed Initial 
Condition technique (|Tormen et al.|1997||Bonafede et al.|20TT]) 



and considering various recipes for the intra-cluster medium 
physics (Fabj an et al.|201 lj Killedar et al.|2012 l. In this work, 
as in R12, we consider only the CSF (cooling-star formation- 
feedback) runs. The finest mass resolutions for the dark matter 
and gas particles are, respectively, m^„, = 8.47 x 10 8 /! _l M Q and 
m-gas = 1.53 x 10 8 !i~ 1 Mq. The Plummer equivalent force used 
to compute the gravitational interaction in the high resolution 
regions is kept fixed to e = 5 h~ l kpc in comoving units from 
redshift z = 60 to z = 2. Thereafter, the same value is kept con- 
stant in physical units. The minimum SPH smoothing length 
allowed is 0.5 x e. 

We assume a flat ACDM model with cosmological param- 



eters in agreement with WMAP7 constraints (Komatsu et al. 



201 1 1: Qm = 0.24 for the matter density parameter, J7a = 0.76 
for the dark energy, Oj ar = 0.04 for the baryon density, Ho = 72 
km s -1 Mpc -1 for the Hubble constant at redshift zero, n s = 0.96 
for the primordial spectral index, and erg = 0.8 for the normal- 
ization of the power spectrum on a scale of 8 h Mpc. 

The hydro-resimulation are performed using the TreePM - 
SPH GADGET-3 code an algorithmically improved version of 
the original code GADGET-2 ( Springel et al. 2005 ). The physics 
implemented in the resimulations analyzed in this paper in- 
cludes a model for heating and cooling from an evolving but 



spatially-uniform UV background (Wiersma et al.|2009i. The 



star formation is modeled using the recipe of |Springel & Hern- 



|quist| ( |2003] >: when the gas particle reaches a density above a 
certain threshold it is treated as a multi-phase particle with the 
hot and cold phases in pressure equilibrium. The last one serves 
as reservoir of star formation. The intracluster medium is en- 
riched by various metals expelled by different stellar popula- 



tions (Tornatore et al. 2007). The galactic winds mimicking 
supernovae explosions have a velocity of 500 km s" 1 . 

The simulation outputs considered here refer to redshift z = 
0.25. In this way, all objects have their /?50(Jj within the field 
of view. The total mass enclosed in /?200 (^500) spans from 4 to 
15 (from 2 to 11) times 10 14 /t'Mq. Other works in literature 
([Poole et al.||20D6l |Jeltema et al][2005l [20081 |Maughan et aL 



2008 ) have shown how the hierarchical formation of structure 
implies a logical change of the distribution of the morpholog- 
ical parameters: a set of high-redshift clusters will be an aver- 
age more disturbed. We remind that, in this paper, we aim at 
testing the efficiency of the investigated parameters rather than 
furnishing precise statistical indications for the distributions of 
the parameters over time. Selecting a medium redshift as 0.25 
is a suitable choice for our investigation. 

2.2. X-ray Catalogue 

To facilitate the post-processing on the X-ray images we ex- 
cluded a priori some particles before creating the event files. 
Namely, we omitted the coldest and densest particles satisfy- 
ing the following condition: T p < 3 x 10 6 /?g a s 5 p- These par- 
ticles exist in the simulation as consequence of the overcool- 
ing phenomenon that typically affects hydrodynamical simula- 
tions with cooling and star formation (see discussion in |Bor-| 
gani & Kravtsov||201 1\ . The normalization value of this re- 
lation was set empirically by looking at the density and tem- 
perature of all particles of each of the twenty clusters. The rela- 
tion was defined such to separate the over-cooled particles from 
the thermalized-cluster particles. The exponent of the density, 
0.25, is linked to the politropic index of the simulations (see 

4 Radius of the sphere whose density is 500 times the critical density at the 
cluster redshift. 
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Appendix A of R12 for details). 

The event files are created thanks to the X-ray MAp Sim- 
ulator (X-M AS) code (presented in |Gardini et "aT|2004| |Rasia| 
et al. 2008 ) to reproduce Chandra-]ike observations. We con- 
sider the ancillary response function (ARF) and redistribution 
matrix function (RMF) of the aim-point of ACIS-S3 through- 
out the field of view of size equal to 16 arcmin (equivalent to 
2561 h~ l kpc accordingly to the assumed cosmology and con- 
sidered redshift, z = 0.25). For the galactic absorption, we adopt 
a WABS model fixing the column density at Nh = 5 X 10~ 20 
cm" 2 . The X-ray emission, instead, has been reproduced us- 



original image: 



F = 



£(/) 



(1) 



ing a MEKAL model in the XSPEC package (Arnaud 



1996) 



and assuming a constant metallicity equal to 0.3 times the so- 
lar metallicity measured by Anders & Grevesse ( 1989 (. The 
exposure time was set to 100 ksec. 

The imaging analysis conducted in this paper, uses packages 
developed within CIAO (Fruscioi ie et al.| |2006) and IDL sub- 
routines. Contrarily to what done in R12, in this work, we do 
not exclude any region or substructure or merging clump within 
7?5oo- The goal of the previous work was, indeed, to measure 
the X-ray hydrostatic mass and therefore we were interested to 
characterize the regularity of the portion of the gas used to com- 
pute the total mass. This paper, instead, aims at characterizing 
the overall X-ray morphology. 

3. MORPHOLOGICAL ESTIMATORS 

In this Section we introduce the morphological parameters 
investigated presenting the original definition as well as their 
modified version adopted in this paper. The first three param- 
eters (Asymmetry, Fluctuation and X-ray Concentration) are 
strongly connected to the definition of the 'CAS' parameters 
provided by Conselice (2003). The original work by Conselice 
was dedicated to optical images. Specifically, the Author was 
looking for a classification method allowing the connection be- 
tween the structural appearance of the stellar light and the galac- 
tic formation history. The acronym 'CAS' stands for Concen- 
tration of light, Asymmetry, and Spatial frequency dumpiness 
parameter. Very similar parameters have been introduced in 
relation to X-ray images throughout the years with modified 
nomenclature: X-ray surface brightness concentration, asym- 
metry, and fluctuation, respectively (e.g. Cassano et al.||2010 
|Okabe et al.|20T0l|Zhang et al.|2010| i. 

Furthermore, in addition to the hardness ratio parameter, we 
test other three indicators (the centroid shift, the third-order 
power ratio, and its maximum). These have been extensively 
used in X-ray works on clusters of galaxies (e. g. |Mohr et al.| 
19951 [Buote & Tsai||19951 |Buote & Xu||19971 |Jeltema et all 
2005] [Bauer et al.|2005||Maughan et al.|2008||VentimigIia et al. 
2008 Jeltema et al.|2008NAndersson et al.|20091|Hallman et al. 
2010 Icassano et al.|2010||Hallman & Jeltema|201 l| |Bohringer 
et al.|2010[|Weifimann et al.|2013| R12). Specifically, two re- 



cent works, have also addressed the issues of Poisson noise and 
X-ray background in evaluating these parameters and offered 
strategies to adopt in real observations (Bohrin ger et al.|2010| 
IWeiBmann et al.|2013| >. 

Fluctuation parameter, F. The fluctuation parameter em- 
phasizes the presence of peaks (valleys) of high (low) X-ray 
flux over the image. Various similar definitions have been his- 
torically used. We prefer an expression similar to the asym- 
metry parameter (see below): F is obtained by subtracting a 
smoothed image, B, from the original one, normalizing by the 



Conselice (2003 ) used a similar parameter, called S, for High- 
Spatial frequency dumpiness that, however, considered only 
the peaks of the surface brightness. Indeed, in the computa- 
tion of the sum in Eq. [T|the absolute value was left out and all 
negative differences of (I-B) were discarded. The parameter 
was computed without the absolute value also by the Locuss 
collaboration ( |Okabe et aT1|2010[ |Zhang et aT71|2010| > but the 
negative deviations were kept in the sum with their sign, in a 
way that positive enhancements were eliminated by negative 
depressions. Their smoothed map was obtained by assuming 
a Gaussian filter of FWHM equal to 2' (equivalent to 400 kpc 
at z=0.2). Accordingly to this assumptions the fluctuation pa- 
rameters varied from to 0.14 and 0.05 was considered the 
boundary dividing regular and disturbed objects. 

We derived three smoothed maps by varying the Gaussian- 
filter FWHM from 2' to 16", and to 8" to test the effect of the 
smoothing on the results. The associated fluctuation parameters 
are, respectively, called F\, Fz, and Ft,. 

Asymmetry parameter, A. One of the most powerful indica- 
tor of the morphology of a galaxy is the asymmetry parameter 
(e.g. |Schade et al.||1995] l obtained by summing the residuals 
between an image, /, and its 180° rotated counterpart, R, nor- 
malized by the original image: 

A ~ ECO () 

In this fashion, all asymmetries non perfectly aligned with the 
diagonals are emphasized. 

The Locuss collaboration applied this definition to soft X-ray 
images ob tained by XMM-Newton ( |Zhang et aT|2010| |Okabe| 
|et al.|2010| ). The images were extracted in the soft band ([0.7- 
2.] keV), corrected for the flat field, point-source-subtracted, 
and binned by 4" x 4". These works pointed out that the clus- 
ters showing an high level of asymmetry are presenting as well 
a discrepancy between the X-ray centroid and weak-lensing 
center. This last feature clearly denotes a recent dynamical ac- 
tivity. The Authors identified as disturbed clusters all objects 
having A > 1.1. The full range of the asymmetry parameter for 
their sample spanned from 0.07 to 1.5. 

In this work we consider three different asymmetry parame- 
ters. Namely, the image R is rotated by 180°, A Iol , ii) flipped 
along the x-axis, A x , and Hi) flipped along the y- axis, A y , with 
both axis passing through the X-ray centroid. Finally, per each 
cluster, we also consider the maximum among the three asym- 
metry parameters (A max ). With this procedure, we avoid any 
possible underestimate of non-asymmetry present along a fa- 
vorite direction. 



Light Concentration, c. Analogously to the definition of 
'light concentration' by Conselice] ( |2003[ l (see also Bershady 
et al.|2000| a specific parameter measuring the concentration of 
the X-ray -emission has been defined by Cassano et al. ( 2010) >. 
This is as: 

S(r < lOOkpc) 

C [kpc] = „ cr ,r.7___ ■ y^i 



S(< 500k P c) 

A similar definition was used by Sant os et al.| ( |2008"| ) as an in- 
dicator for the presence of cooling core systems at high red- 
shift: c = S(r < 40kpc)/S(r < 400kpc). The radii of 40 and 400 
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kpc were chosen to maximize the separation of concentration 
values between cool-core (cc) and non cool-core (ncc) clus- 
ters (see also |Semler et al. :2()l 2). However, simulations with 
radiative cooling and not powerful feedback such as that pro- 
vided by active galactic nuclei are still unable to reproduce the 
statistic of cc versus ncc systems (see the review by Bor gani &| 
Kravtsov 2011, regarding the over-cooling problem affecting 



with respect to the center. The parameter is defined as follow: 



simulations). For this reason, we select a larger region to repre- 
sent the core. Hall man & Jeltema| ( |201 1) noted the limitations 
of considering apertures at fixed physical radii in the study of 
large samples of clusters at different redshifts. They proposed 
the usage of apertures with radii expressed in units of R500 as 
done in Kay et al. ( 2007| l. Taking into account this comment, 
in addition to the definition presented in Eq. [3] we test another 
definition of the concentration: 

_ S(r < 0.2 x R 500 ) 

C [R500] - Cf , p S • W 

i>{r < K500J 

Cassano et al. ( 20i0| ) studying 32 clusters observed with Chan- 



dra (with redshift between 0.2 and 0.4 and X-ray luminosity 
above 10 44 erg s" 1 ) proposed Ck pc = 0.2 (median value of their 
distribution spanning from ~ 0.05 to ~ 0.7) as the dividing line 
from regular and disturbed systems. 

Hardness ratio indicators, H. In the hierarchical structure 
formation scenario larger system form through accretion of smaller 
structures. In this picture, a merger only sporadically involves 
two objects with the same mass and, most likely, happens be- 
tween two or more systems presenting a large spread in mass 
that can be translated into a spread in temperature. A good 
variety of movies based on simulations clearly shows how the 
thermalization process of the intra-cluster gas is continuously 
challenged by the incursion of smaller and colder blobs |^| The 
multi-temperature phase of the ICM of a merging object can be 
enhanced by taking the ratio between images obtained in the 
soft band (such as [0.3-1.5] keV where the cold gas is strongly 
emitting) and in the hard band (such as [1.5-7.5] keV where the 



hot gas has more power). This technique, recently used by Gitti 
|et al.| ( |201 1) to identify cold filaments in Hydra A cluster, has 
a long history of usage to identify temperature structures (e.g. 
|Allen & Fabian|1997[|Ettori et al.|1998) . 
In R12, we considered two hardness ratio 



Z(H-S) 
S(S) 



and H 2 =T,(H/S) 



(5) 



where H and S are the images in the aforementioned hard and 
soft energy bands, respectively. In this paper, we focus only on 
the first definition that has been proved to be more effective and 
we also consider another ratio defined with S and H extracted in 
the soft energy band [0.7-2] keV and in the hard energy band [2- 
7] keV, respectively. Before computing the ratios, we applied to 
both images a Gaussian filter of FWHM equivalent to 2', 16", 
and 8". 

Centroid shift, w. This measure, as the name suggests, as- 
sesses how much the centroid of the X-ray surface brightness 
moves when the aperture used to compute it decrease from a 
certain R max to smaller radii. The presence of X-ray bright 
clumps can produce significant changes on the X-ray centroid 
unless they are distributed in a perfectly symmetric geometry 



5 See, for example, the movies at the following link http://www.mpa- 
garching.mpg.de/galform/data_vis created by Klaus Dolag 
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W = 



£(A- < A >) 2 



(6) 



where N is the total number of apertures considered and A, 
is the separation of the centroids computed within R mdx and 
within the i th aperture. This parameter has been extensively 
used in literature (e.g. |Q'Hara et al.||2006| |Poole et aL[|2 006; 
Maughan et al.||2008} [V entimigri a et al.||2008[ |Jeltema et aL| 
2008||Bohringer et al'|2010||Weifimann et al.|2013| > since it is 
an useful characterization of the dynamical state of the cluster 
(Mo hr et al.|1993| l. The robustness holds also if the parameter 



is slightly differently defined ( |Thomas et al.|[1998t Kay et al. 
2007 ). We compute Eq.[6]by varying the apertures from 0.15x 
^max> to ,R m ax with a step of 0.05 x R max and R max varying from 

30% to 1 #500. 

The most recently chosen boundary between regular and dis- 
turbed objects present in literature is around 0.01. Specifically, 
w = 0.012 for Cas sano et al.| ( |20T0) l who studied 45 luminous 
clusters observed with Chandra, w = 0.01 for WeiBm ann et al.| 
pOH]) for 80 cluste rs observed with XMM-Newton, 0.02 for 
O'Hara et al. (2006 ) who analyzed 45 clusters from the ROSAT- 
PCPS images (thus with less spatial resolution than Chandra). 
The ranges reported in literature span from ~ 0.001 to ~ 0.15 
([Poole et al.|2006||Ventimigha et al.p008l|Jeltema et al.|2 008; 
|Cassano et al. 1 ,2010; WeiB mann et al.||2013| ) with some small 
departures on the limits depending on the R max chosen. 



Power ratios, P3 /Pq. The power ratios, introduced by Buote 
|& Tsai ( |1995| l, mimic a multiple decomposition of the two- 
dimensional projected mass distribution inside a certain aper- 
ture, R dp . Instead of the mass, however, they are applied to the 
X-ray surface brightness images, S. 

The generic m-order power ratio (m > 0) is defined as P m /Po 
with 

P ™= 2m 2R2m ( a m +b l) and Po = a Q \n(R. dp ) (7) 

where ao is the total intensity within the aperture radius. The 
generic moments a m and b„, are expressed in polar coordinates, 
R and 0, and given by 

a m (r)= / S(x')R'cos(m(t>')d 2 x', (8) 

jR'<R m 



and 



b m (r)= I S(x')R'sin(m(j)')d 2 x'. (9) 



The quadrupole power f_ refers to the ellipticity of the clus- 
ters, P3 informs about bimodal distribution, P4 is similar to Pj 
but more sensible to smaller scales (indeed, Pi is strongly cor- 
related with P\). The third order power ratio is the most suitable 
to identify asymmetries or presence of substructures. This is the 
parameter investigated in this paper. In literature, various aper- 
ture radii have been used both at a precise physical scale (either 
0.5 Mpc or 1 Mpc) or at #500- In this paper, we will test all these 
possibilities as well as the new parameter described in the next 
paragraph. In most recent literature , the decimal logarithm of 
P 3 /granges from ~ -8.7 to ~ -5.7 (jjeltema et al |20 05 2008 ; 
|Cassano et al.||2010} [WeiBmann et al. 2013} again depending 
on the maximum radius used for the calculation of the power 
ratios. The median values are typically logio(i_/./_) = _ 7.2 and 
the dividing limit between regular and disturbed system is set 
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around -6.7 and -7. WeiBmann et al. ( 201 3[ > proposed other 
two boundaries to select exclusively regular and disturbed ob- 
jects: Pi/Po< 10" 8 andP 3 /P > 5 x 10" 7 , respectively. 



1 x R 



Maxi mum of the third pow er ratio, P 3 /P 

paper. 



In a recent 



WeiBmann et al. (2013 i discussed the limits of the third 
order power ratio computed within a fixed radius (either in phys- 
ical units or in units of /?5oo). Their claim is based on the con- 
sideration that the P3/P0 of the Bullet cluster evaluated at R500 
has a value typical of a regular system (~ 10~ 7 ). This aston- 
ishing result is explained by the fact that the influence of the 
'bullet', located in a more central zone, is smeared out when 
considering the entire cluster volume. The evaluation of the 
third order power ratio at a fixed radius limits the power to dis- 
criminate between 'regular' and 'disturbed' objects. To fixed 
this problem, WeiBman n et al.| ( [2013) > proposed to consider the 
maximum of the third order power ratios estimated at various 
radius from 0.3 to 1 ^500 with step equal to 0.1 x ^500- We, 
therefore, compute all eight power ratios, and per each cluster 
consider the maximum among them. 

4. RESULTS 

Unless otherwise specified, all morphological parameters are 
measured on the soft X-ray image binned 2" x 2" at various 
ratios of/?goo: 0.3,0.5,0.8, andl. The efficacy is proved against 
the visual classification established after a careful process. 

Visual classification. In order to closely reproduce the obser- 
vational approach we divide the sample into three classes: reg- 
ular, disturbed, and 'mix'. To avoid influence of the observer- 
bias, three different people looked at and rated the images. These 
were shown three times each and appeared in random order. For 
every image, we average the nine grades. These span from 
(given to a spherical system without substructures) to 3 (very 
disturbed X-ray emission or clear evidence of large merging 
substructures). Strongly elliptical clusters or spherical clus- 
ters with minor substructures still in the outskirts (that have 
not crossed R5QQ) and not interfering with the X-ray emission 
of the main system resulted typically rated from 1 to 2. All 
images with an average grade below or equal 1.5 are consid- 
ered of regular systems, while those with a grade above 2.5 are 
proper of disturbed objects. The two subsamples are shown in 
Fig. [8] and Fig. [9] of the Appendix. All other objects fall into 



a class we call 'mix' (Fig. 10 1. As evident from the images, 
our regular clusters are part of the 'single', 'elliptical' (even 
if not strongly elliptical), and 'galaxy' categories proposed by 
|Jones & Forman| (1992) while the disturbed system could be 
represented by their 'complex', 'double', 'off-center' classes. 
In next Section, we will study the efficiency and improvement 
of combining more estimators together and will calculate some 
statistical parameters as purity and completeness. For this rea- 
son, we will need to consider the whole sample. Consequently, 
we proceed to subdivide the 'mix' class into semi-regular and 
semi-disturbed. The limit of the averaged grade chosen to di- 
vide the two extended classes is set equal to 2. 

Asymmetry and fluctuation parameters We investigate the 
efficiency of the asymmetry and fluctuation parameters in dis- 
tinguishing among regular and disturbed clusters both individ- 
ually and combining them together (Okab e~et al.|2010" Zhang 
|et al.|2010) . 
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FIG. 1. — Top panel: distribution of the regular (red) and disturbed (blue) 
clusters as a function of the fluctuation parameter, F. The smoothed image 
used in Eq. [T] is derived with a Gaussian smoothing of FWHM equal to 2' 
(F[, upper panels), 16" pixel (F^, middle panels), and 8" pixel (F3, bottom 
panels). The sum is computed over a circular region of radius specified at 
the top of the columns. Bottom panel: variation of the fluctuation parameter 
statistic (vertical values show the range between and 172 °f Table 1). The 
numbers listed within the Figure show the fraction of R500 used in the sum of 
Eq.[T] 



Per each of the three smoothed maps, we evaluate F at four 
different radii, and then compare all results with the visual clas- 
sification. In the left panel of Fig. [T] we show the histograms 
of F\,Fi, F3 for regular (in red) and disturbed (in blue) clusters 
(for clarity, we omit the 'mix' class). The distributions of the 
whole sample, regular sub-sample and disturbed one are sta- 
tistically characterized by the median, the first and third quar- 
tiles, reported in Table[TJand plotted in the right panel of Fig. [TJ 
(for the following parameters we will omit the figure since the 
information is already contained in the tables). We, further, ap- 
ply a Kuiper-two-sample-test statistics to determine the prob- 
ability that the two distributions are statistically different. For 
the numbers in this work (20 regular clusters and 13 disturbed 
ones) a probability p below 0.05 implies a net distinction be- 
tween the two distributions. 

We conclude that the fluctuation parameter is a good indi- 
cator of the X-ray morphology, in particular, when the mea- 
sures are done at Rsoa and the largest Gaussian filter is applied 
{Fx performs better than F^ in all cases apart in the real cen- 
tral zone). Specifically, the limit of F\(?A /?5oo) = 0.9 effectively 
divided the two morphological samples. Measuring the param- 
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eter at the center of the cluster (0.3 x R500) is a good strategy to 
select regular objects: when Fi < 0.6 (F 3 < 0.5) 84% (82%) of 
the objects selected are regular and the medians of the distri- 
butions are significantly apart. 
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FIG. 2. — Top panel: distribution of the regular (red) and disturbed (blue) 
clusters as a function of the asymmetry parameters: A rot (top panels), A x (cen- 
tral panels), and A max (bottom panels). The sum is computed over a circular 
region of radius specified at the top of the columns. Lower panel: distributions 
of fluctuation (Fi) and asymmetry (A ro t) parameters for all the clusters. The 
values of horizontal and vertical lines are listed in TablefS] 

We assess the performances of the asymmetry parameters 
following the same procedure. The results for the entire analy- 
sis are reported in the bottom panel of the table (Table [TJ. The 
distribution of A for regular and disturbed objects are shown 
in Fig. [2] only in the cases of rotation (top panels), of flipping 
around the x— axes (central panels), and of the maximum of all 
asymmetry parameters. The case A y is extremely close A x and, 
therefore, we chose to avoid its listing in Table[JJ The asymme- 
try parameter distinguishes effectively between the two classes 
of clusters, in particular when Eq. [2] is derived including large 
radii. Selecting clusters with A ro , < 1.1 at R500 only one dis- 
turbed cluster (5% of the sample) contaminates the sample of 
otherwise regular objects. The two distributions are generally 
broad and often present two major peaks (this is especially true 
for the regular systems). The most notably advantage of using 
the maximum of the three asymmetry parameters is visible in 
the bottom-right panel where the two distributions show clearly 
separated peaks with small dispersion (see also Table [TJ. Ex- 
cluding this case, A max does not significantly improve the clas- 



sification. For this reason, in the rest of the paper we will focus 
only on A mt computed at R500, the version most used in litera- 
ture. 

Selecting the clusters by combining the constraints on both 
the asymmetry and the fluctuation parameters, A rot (at Rsoo)< 
1.15 and Fi(at R500) < 1-0, quite robustly gives regular objects 
(only one clear contaminanrl) while objects with A > 1.15 and 
F > 1 .0 are mostly disturbed with no contaminants (see right 
panel of Fig. |2j. We conclude by noticing that the asymmetry 
parameter covers the same range of values as those derived by 
the Locuss collaboration (Zhang et al. 2010| Okabe et al. 201Q) >. 
On the other hand, we cannot compare our and their fluctua- 
tion parameters since F is very sensitive to the area over which 
the sum of Eq. 1 extends and none of these papers specifies 
the maximum external radius used. For our simulated sample, 
adopting their formula we obtained F values from to 0.03 
when we integrate over R 5 qq and from 0.03 to 0.14 when we 
stop at half of that radius. 

X-ray SB concentration The two definitions of X-ray surface- 
brightness concentrations {c[kpc\ an d c tR 5 oo]) are quite different 
since the R^q values of our clusters range from ~ 1 Mpc to 
~ 1.5Mpc (see Table 1 of R12 where the radii are reported in 
/r'kpc). The results related to this parameter are reported in 
Tableland shown in Fig. [3] When the measurements are per- 
formed using the physical radii the two distributions have very 
small dispersions and two clearly distinct peaks that, however, 
are as close as to produce an overlap of the tails of the distribu- 
tions (left panel). The histograms significantly shift apart when 
the two radii used to compute the surface brightness are in units 
of /?5oo (right panel). Both approaches are well suited to create 
subsamples of objects. In particular, in the first case, none of the 
objects with c^ pc < 0.15 (c^ > 0.2) are regular {disturbed), in 
the second case a similar situation is presented for Cr 500 < 0.25 
(c« soo > 0.3). Even so, there is a significant difference between 
the numbers of objects that are located in between these bound- 
aries: 4 regular (20% of the regular sample) and 7 disturbed 
objects (~ 50% of the disturbed sample) have Ck pc between 
0.15 and 0.2, while only 2 regular (10%) and 2 disturbed sys- 
tems (15%) have cr 500 between 0.25 and 0.3. From Table|2j we 
notice that our overall distribution has the same median of the 



Chandra sample studied in Cassano et al. (2010). 



Hardness ratio indicators We find almost no difference be- 
tween the distributions of regular and disturbed clusters classi- 
fied on the basis of the two hardness ratio parameters, H, for 
all the smoothing lengths applied and for all radial range con- 
sidered. The p value of the Kuiper-two-samples test associated 
to the two H distributions of regular and disturbed objects is 
always above 0.055 meaning that the this parameter is not effec- 
tive in distinguishing the two morphological classes. In Fig. |4] 
we present two examples of H (with S extracted in the [0.3- 
1.5] keV band and H in the [1.5-7.5] keV band) computed after 
smoothing the images with a Gaussian of FWHM equal to 16" 
and extending the sum to the smallest (R < 0.3 x ^500) and to 
the largest (R < R500) circular region. In these cases, p = 0.187 
and p = 0.057, respectively. All other histograms (including 
those obtained from a different definition of S and H) look very 
similar: the disturbed-clusters are characterized by two peaks 

6 We call contaminant a cluster non part of the 'mix' class and lying in the 
parameter space typical of the other morphological class. 
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Table 1 

Fluctuation and asymmetry parameter distributions (median, first and third quartiles) computed at various 
fraction of r500 for the regular sample, all the clusters, and the disturbed sample. last column refers to the 

kuiper-two-sample-test probability. 



par 




regular 






all 






disturbed 










q lst 


med 




q Ut 


med 


q 


3rd 


q lst 


med 


q 


3rd 




p 


Fx (R500) 


0.81 


0.87 


0.93 


0.80 


0.90 


1 


01 


1.01 


1.02 


1 


06 





008 


Fi (R500) 


0.61 


0.72 


0.83 


0.65 


0.76 





91 


0.90 


0.93 





97 





001 


F3 (^500) 


0.33 


0.38 


0.46 


0.63 


0.73 





87 


0.48 


0.59 





64 





001 


F x (0.8 x R 500 ) 


0.77 


0.82 


0.88 


0.76 


0.84 





95 


0.94 


0.96 


1 


02 





012 


F 2 (0.8 x R 500 ) 


0.58 


0.67 


0.78 


0.60 


0.72 





85 


0.84 


0.88 





92 





012 


F 3 (0.8 x R 5QQ ) 


0.55 


0.64 


0.75 


0.58 


0.69 





82 


0.80 


0.85 





88 





012 


F x (0.5 x R 500 ) 


0.64 


0.68 


0.71 


0.64 


0.68 





76 


0.74 


0.81 





87 





012 


F 2 (0.5 x R 500 ) 


0.46 


0.53 


0.62 


0.48 


0.60 





70 


0.66 


0.75 





81 





033 


F 3 (0.5 x R 500 ) 


0.44 


0.51 


0.59 


0.46 


0.57 





67 






u 


78 





033 


Fx (0.3 x R 500 ) 


0.52 


0.54 


0.57 


0.52 


0.56 





62 


0.60 


0.68 





71 





033 


F 2 (0.3 x R 500 ) 


0.34 


0.40 


0.48 


0.36 


0.46 





54 


0.50 


0.61 





66 





004 


F 3 (03xR 500 ) 


0.33 


0.38 


0.46 


0.34 


0.45 





52 


0.48 


0.59 





64 





001 


A r {R50Q) 


0.84 


0.93 


1.08 


0.87 


1.06 


1 


17 


1.29 


1.30 


1 


34 


2 


E-4 


A x (R500) 


0.83 


0.93 


1.08 


0.87 


1.07 


1 


18 


1.19 


1.30 


1 


37 


6 


E-4 


A r (0.8 x R 5Q0 ) 


0.80 


0.89 


1.06 


0.83 


1.02 


1 


10 


1.23 


1.25 


1 


29 


2 


E-4 


A x (0.8 x R 5m ) 


0.80 


0.89 


1.05 


0.83 


1.03 


1 


14 


1.14 


1.27 


1 


31 





004 


A r (0.5 x R 5Q0 ) 


0.69 


0.77 


0.92 


0.74 


0.89 





96 


1.05 


1.08 


1 


12 





004 


A x (0.5 x R 5m ) 


0.68 


0.77 


0.90 


0.73 


0.88 





99 


0.99 


1.10 


1 


12 





020 


A r (0.3 x fl 5 oo) 


0.57 


0.64 


0.76 


0.59 


0.69 





83 


0.85 


0.88 





91 





008 


A x (0.3 x fl 5 oo) 


0.55 


0.64 


0.75 


0.59 


0.72 





84 


0.83 


0.88 





92 





021 


A m ax (^500) 


0.85 


0.94 


1.10 


0.89 


1.09 


1 


18 


1.29 


1.34 


1 


37 


8 


E-5 


A max (0.8 X #500) 


0.81 


0.90 


1.06 


0.85 


1.05 


1 


14 


1.25 


1.28 


1 


31 


3 


E-4 


A max (0.5 x R 5Q0 ) 


0.70 


0.78 


0.93 


0.75 


0.91 





99 


1.06 


1.11 


1 


14 
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A max (0.3 x fl 5 oo) 


0.57 


0.66 


0.78 


0.60 


0.74 





85 


0.85 


0.91 





95 
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Table 2 

x-ray surface brightness concentration, centroid shift, and third order power ratio parameter distributions 
(median, first and third quartiles) for the regular sample, all the clusters, and the disturbed sample. last 

column refers to the kuiper-two-sample-test probability. 



par 




regular 






all 






disturbed 








H 


nicd 


H 




mcd 




H 


nicd 


a 3ld 

H 




SB(0.2 x R 5O0 )/SB(R 500 ) 


0.35 


0.42 


0.48 


0.27 


0.37 


0.44 


0.19 


0.21 


0.23 


7E-4 


SB(100kpc)/Sfi(500kpc) 


0.20 


0.24 


0.25 


0.16 


0.20 


0.25 


0.14 


0.15 


0.16 


0.004 


l°gio( w ) 
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9 HQ 

-z.uy 


1 79 
- 1 . /Z 


9 HQ 

-z.uy 


1 78 
-1. /o 


1 AQ 


1 

-1 .Do 


1 9H 
- 1 .ZW 


i in 

-1 . 1U 


H H9H 


^3 1^500 ) 


7 *\A 


7 91 

- / .Z 1 


6 QzL 


7 1 ^ 
- / . 1 j 


6 7S 
-u. / J 


-U. 1 J 


S 7^ 
-J. / J 


-J.JJ 


A Q7 


9 F S 


P 3 (0.5 x R 5Q0 ) 


-7.32 


-6.95 


-6.48 


-7.06 


-6.52 


-5.84 


-5.79 


-5.59 


-5.34 


0.001 


P 3 (0A x R 500 ) 


-7.01 


-6.64 


-6.38 


-6.77 


-6.45 


-5.93 


-6.34 


-5.63 


-5.47 


0.127 


P 3 (0.3x/?500) 


-6.86 


-6.59 


-6.38 


-6.84 


-6.48 


-6.05 


-6.61 


-6.29 


-5.77 


0.168 


P 3 (500 kpc) 


-7.51 


-7.09 


-6.53 


-7.29 


-6.53 


-5.84 


-5.72 


-5.32 


-5.21 


4E-5 


P 3 (1000 kpc) 


-7.65 


-7.18 


-7.03 


-7.11 


-6.75 


-6.17 


-5.82 


-5.26 


-5.00 


8E-5 


^3, max 


-6.65 


-6.45 


-6.26 


-6.41 


-6.17 


-5.60 


-5.47 


-5.13 


-4.76 


7E-5 
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FIG. 3. — Distribution of the regular (red) and disturbed (blue) clusters as 
a function of the X-ray surface brightness concentration ew on the left and 
cr 500 on the right. 
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FIG. 4. — Distribution of the regular (red) and disturbed (blue) clusters 
as a function of the fluctuation parameter, H. The smoothing image used in 
Eq.|5]assume a Gaussian filter of FWHM equal 16". Only results referring 
to circular regions of radius 0.3xK5oo(left panel) and i?5oo( rl ght panel) are 
shown. 



one of which is always at low H values and, therefore, within 
the range typical of regular objects. Due to the no-success of 
the hardness ratio indicator, we are not reporting the median 
and quartiles of the distributions and we will not consider it any 
further. 

Centroid-shift The histograms of the centroid shift for the 
regular and disturbed clusters present two well separated peaks 
in case of R mdx = /?soo(Fig. [5]). None of the objects having 
w < 0.015 (w > 0.04) is disturbed (regular). Nonetheless, 
between 0.015 and 0.04 the distributions overlap significantly. 
In this range, we find 30% of the total distribution of the dis- 



turbed systems and 35% of the regular ones. Therefore, even 
if the centroid shift is efficient in selecting systems that are part 
of only one category (the median are quite far apart) it also 
presents significant contamination. The ranges of the centroid- 
shift obtained (from 0.002 to 0.12 ) agree well with the val- 
ues from sample observed with C/!ant/ra,XMM-Newton, and 
ROSAT. 

Subsequently, we allow /? max to vary from R500 to 0.3 x Rsoq 
with a step of 0.1 x R500, We found that the Kuiper-two-sample 
probability increases and soon becomes unacceptable when the 
maximum radius decreases below 0.7 x R 5{m . 
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FIG. 5. — Distribution of the regular (red) and disturbed (blue) clusters as a 
function of the decimal logarithm of the centroid-shift parameter, log 10 (w). 
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FIG. 6. — Distribution of the regular (red) and disturbed (blue) clusters as a 
function of the third-order power ratio computed within R500 (right panel) and 
maximum of third order power ratio (left panel). 



Third-order power ratio and its maximum We compute 
the third-order power ratio in eight different apertures: from 
0.3 x /?5oo to /?5oo with step of 0.1 x R500. For each of them 
we characterize the distributions of regular and disturbed sys- 
tems and report the values of medians and quartiles in Table [2] 
The performances of the power ratio quite depend on the radius 
considered. The best results are obtained for a large area (with 
R > 0.5 x /?5oo). Restricting the analysis to the central region, 
instead, impacts negatively on the efficiency of the method. In- 
deed, the p value jumps from 0.001 (amply satisfying the re- 
quirement of p < 0.05) at 0.5 x ^500 to 0. 127 (meaning that the 
distributions of regular and disturbed systems are not clearly 
distinguishable) at lower radius. On the left panel of Fig. [6] we 
show the result of the third order power ratio evaluated at /?5oo- 
The two distributions are significantly separated and no regular 
(disturbed) objects are found with P3/P0 above 2 x 10~ 7 (be- 
low 5 x 10~ 7 ). Similar behavior is found when the power ratio 
is measured within a fixed physical scale such as 500 or 1000 
kpc (see Table. [2]). Also in this case, we observed that measure- 
ments conducted over a larger area are more effective. 

Subsequently, per each cluster we consider the maximum of 
the third-order power ratio profile (i.e. the maximum of the 
eight values computed) as suggested in WeiBmann et al. (2013 1. 
This case is shown in the right panel of Fig.[6]and reported at the 
bottom of Table|2] Per construction, all the numbers are shifted 
towards higher values causing an approach of the two medians. 
The two distributions are narrower and present higher peaks. 
No regular objects are found with P3. max > 3 x 10~ 6 while only 
two disturbed systems are found below this limits (implying a 
9% of contamination). 

5. PAIRS OF PARAMETERS 

In this Section, we study the combination of different pairs 
of indicators and test their robustness as we did for the asym- 
metry and the fluctuation parameters (Fig [2}. All planes of the 
pairs studied are divided in four quadrants by the limits, L, listed 
in Table [4] The parameters used are: the asymmetry parameter 
A rot „ the decimal logarithm of the third order power ratio evalu- 
ated within /?5oo, the X-ray concentration cr 500 , and the decimal 
logarithm of the centroid-shift. The other parameters (asym- 
metry computed by flipping the image around either a vertical 
or an horizontal axes, hardness ratios, third order power ratios 



Table 3 

Coefficients of the linear fit to the points of various 

INDICATOR PAIRS (I s1 AND 2 nd COLUMNS) AND THE PEARSON 
CORRELATION COEFFICIENT (3 rd COLUMN). 



par 


linear fit 


corr. 


pairs 


A 


B 


coef. 


log 10 (P 3 ),A 


0.15 


2.05 


0.72 


logio(P 3 ),c 


-0.08 


-0.18 


-0.73 


logio(P3),logio(w) 


0.27 


-0.001 


0.66 


c,A 


-1.62 


1.63 


-0.87 


logio(w),A 


0.31 


1.61 


0.66 


logio(w),c 


-0.19, 


0.008 


-0.77 


A mU F 


0.55 


0.33 


0.85 



Table 4 

Per each pair of indicator: values of the limit (L) 

DIVIDING REGULAR AND DISTURBED OBIECTS (VERTICAL AND 
HORIZONTAL IN FlG.[7} AND THE PURITY (P) AND 
COMPLETENESS (C) PARAMETERS. 



par 




L 


regular 


disturbed 


pairs 


ver. 


hor. 


Pr 


cv 


Pd c d 


log 10 (P 3 ),A 


-6.4 


1.15 


0.95 


0.95 


1.0 0.77 


logio(P3),c 


-6.4 


0.32 


1.00 


0.87 


0.9 0.45 


logio(^3),loglo(w) 


-6.4 


-1.4 


0.95 


0.95 


1.0 0.45 


c,A 


0.32 


1.15 


1.00 


0.87 


1.0 0.72 


logio(w),A 


-1.4 


1.15 


0.93 


1.00 


1.0 0.36 


log 10 (w),c 


-1.4 


0.32 


0.97 


0.87 


1.0 0.45 


A ro t,^ 


1.15 


1.00 


0.89 


0.89 


1.0 0.83 



within smaller circular regions, and the maximum of the third 
order power ratio) are not considered here since their perfor- 
mances are either worse or similar to the selected indicators. 
The results are shown in Fig. [7] Per each combination, we also 
compute the linear fit and the Pearson correlation coefficient 
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(Table |3J. All pairs show a significant correlation degree (cor- 
relation coefficient always above 0.66). The pairs that correlate 
the least are those involving the centroid shift that for our sam- 
ple poorly accomplish the result of dividing the two samples 
and present a large scatter of values (see Pig . |5j . The pairs that, 
instead, correlate the most are those combining the asymmetry 
parameter with either the concentration or the fluctuation pa- 
rameters. However, looking at Fig. [2] and Fig. [JJ we can notice 
that several objects (~ 5) clearly identified either as regular or 
disturbed are not located in the quadrant defining their mor- 
phological class. 

To establish the best combination of parameters, we, there- 
fore, calculate the 'purity' and 'completeness' of the regular 
and disturbed classes. To enlarge the statistics, we include the 
'mix' class spitting it into two other classes: semi- regular and 
semi-disturbed (see Appendix). In Fig. [7] these two categories 
are shown by magenta and green points, respectively. We pro- 
ceed then to evaluate the purity and completeness of the two 
extended classes. In the case of the regular plus semi-regular 
systems (20+18 clusters) the two parameters are, respectively, 
defined as: 



QN(regular) 



and Cr Regular) (1Q) 



QN(disturbed+ regular) TN(regular) 

where QN represents the number of objects in the quadrant 
and TN the total number of objects. In similar way, we also 
compute Pj and Q with respect to the disturbed plus semi- 
disturbed class (13+9 objects). The resulting values are re- 
ported in Table [4] together with the limits L chosen per each 
parameter. 

From this analysis, several interesting conclusions emerge. 
Purity of both morphological classes is very high (always above 
95% except two cases) implying that there is very little contam- 
ination when two indicators are combined. The double bound- 
aries, therefore, clearly enhance the power of the estimators in 
distinguishing between the two morphological classes. In most 
case, nevertheless, a good number of objects is lost, especially 
in the case of the disturbed sample. Premise that this result 
might be influenced by the smaller number of disturbed objects 
(the lost of four systems reduces Q to ~80%), we notice that 
completeness is particularly poor whenever the centroid-shift 
parameter is involved (45%, 36%, and 45% for the combina- 
tion with log 1( )(P3/P ), A, c, respectively). Taking into account 
all the parameters investigated the best combinations involve 
the asymmetry parameter and either the third order power ratio 
or the concentration parameter. 

Encouraged by the improvement reached by joining two indi- 
cators, as final step, we defined a new parameter that combine 
a larger number of parameters. Considering the ones studied 
in this Section (par = [logio(Pj,/Po),logio(w),c,A,F]), M is as 
follows: 



M - 



1 



par -L par 



P ar nar \Hpar " L p 



(11) 



par 



Np a r is the number of parameters analyzed (N par = 5 in our 
case). A p ar regulates the sign: it is considered equal to +1 only 
in the case of par = c, otherwise it is fixed equal to -1. In this 
fashion, regular systems should have a positive M. L par and 
m par are, respectively, the limits used in Fig. 7j (and listed in 
Table |4| and the medians of the distribution ofthe parameters 
par extended to the whole sample (Table [Tj and Table [2]). Fi- 
nally, q par is either the first or the third quartile depending if 
the parameter of the specific cluster is smaller or larger than its 



median, m p . In other words, M is the sum of all differences 
of the parameters from their limits normalized by the distance 
from quartiles to the medians. In the left panel of Fig. [7] we are 
showing the histograms of M. Red and blue distributions repre- 
sent regular and disturbed objects while red-magenta and blu- 
green ones include also semi-regular and semi-disturbed sys- 
tems. The separation between the two more restrict morpholog- 
ical classes is clear, and the overlap for the extended classes is 
limited to three objects, all close to the boundary. All objects 
with M > 0.05(< -0.05) are regular {disturbed). As evident 
from the image, the Kuiper-two-sample-test probability, in both 
cases, is extremely low (p = 8 E-5 and p =1 E-8, respectively). 
One of the most remarkable feature is that the majority of the 
objects in both morphological classes are located in regions sig- 
nificantly apart from 0. This is especially true for the regular 
clusters whose distribution's peak is located around 2, meaning 
that most of the clusters have more than one parameter signifi- 
cantly beyond the quartile of its distribution. From comparing 
the distribution plots in Figures 1 through 5 with the quartiles 
values listed in Table 1 and 2, we evince that the M parameter 
is the only one presenting this characteristics. For all the other 
parameters, indeed, the peak of the regular clusters' distribu- 
tion is around to the value of the quartile ( M ~ 1). Since the 
peaks of the two distributions are far apart from the boundary 
and greatly apart from each other, M is the strongest measure of 
the X-ray morphology. We further compute M assuming N par 
equal to 4 (removing the centroid shift) and to 3 (removing also 
the fluctuation parameter). The distributions in both cases are 
very similar to the ones plotted in Fig. [7J (results not shown). 
However, we noticed that the two morphological classes have 
higher peaks when we use all five indicators. The only advan- 
tage we witness by excluding the centroid-shift parameter is 
that the peak of the disturbed class moves farther away towards 
more negative numbers. The last but probably most important 
feature of the M parameter is that it is equivalently strong in dis- 
tinguishing the two morphological classes even when in Eq 11 
we use the median of each parameters instead of the limit L par . 
Also in this context the distribution of the total sample is clearly 
binomial (such as the distribution shown in the right panel of 
Figure 7). This property is not present in any of the individ- 
ual parameter for which none of their median is equivalent to 
the boundary, L par , individuated in this work. Since M efficacy 
does not depend on a a priori knowledge of the boundaries, it 
is the most accessible tool to use in future surveys (provided 
that the objects observed are numerous and not biased towards 
a particular morphological class). 

6. CONCLUSIONS 

Analyzing 60 synthetic Chandra images, we tested several 
morphological parameters commonly used in the literature to 
distinguish between regular and disturbed systems. The mor- 
phology of each image has been rated 9 times (three grades 
assigned by each of the three observers) with a grade from 
to 3, where indicates extremely regular objects and 3 very 
disturbed ones. All systems with an average grade < 1.5 or 
> 2.5 were declared as regular or disturbed, respectively. In 
the final count the two classes have 20 and 13 objects. We used 
this definition to evaluate the capabilities of the estimators to 
identify systems clearly regular or clearly disturbed. Subse- 
quently, we studied various combination of parameters. For 
this goal, we extend the analysis to the whole sample to enlarge 
the statistics and to be able to compute the purity and complete- 
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ness parameter. We proceeded to sub-divide the 'mix' class into 
semi-regular and semi-disturbed objects (see Appendix) fixing 
the diving grade equal to 2. We decided to test the morpho- 
logical parameters against the visual classification because our 
intent is to search for the best parameter that can be as accurate 
as a trained eye. However, a posteriori, we checked also at the 
mass accretion history of our objects and found that our regular 
and semi- regular clusters have accreted on average only 10% 
of their mass at ^200 m the last megayear against the 20% of the 
most disturbed systems. These numbers grow, respectively, to 
30% and 50% with respect to a 1.7 Myr time lapse. Looking at 
the mass accretion history of each individual cluster we found 
that all our regular systems are growing their mass slowly with 
the only exception of CL20. This system experienced an almost 
simultaneous merger by three different structures that, however, 
are still outside the Chandra field of view. As a confirmation all 
the morphological parameters investigated returned values typ- 
ical of a regular object. 

For each estimator, we showed the distribution of the two ex- 
treme morphological classes and reported the median, 1 st and 
3 quartile. We discussed the level of performance of each in- 
dicators by listing the percent of contaminant systems (defined 
as clusters of a particular morphological class lying in the pa- 
rameter range typical of the other class). Finally, we analyze the 
advantages of considering pairs of indicators over a single one. 
Following this logic, we also define a new parameter, M, that 
adds the performances of five indicators. This new parameter is 
the strongest measure to discriminate between morphological 
classes among those tested in this paper. 

Before summarizing our results, let us remind that the val- 
ues of all the parameters computed using our synthetic cata- 
logue agree with those derived from cluster samples based on 
observations apart the X-ray surface-brightness concentration 
parameter. Indeed, while the median value of this parameter is 
coincident with that found in Cas sano et al.| ( |20fD) , our concen- 
trations reaches at maximum the value of 0.4 while those from 
the observational sample extend to c = 0.7. We advance two 
plausible explanations for this difference: z) our clusters are all 
very massive and if) the majority of them lie in a dense envi- 
ronment (R12). The first fact leads to a lower concentration of 
mass (and, thus, light) at the center of the objects with respect 
to smaller groups (Nav arro et al.||T996 1. The second situation 
implies that out systems have a large amount of mass in their 
outskirts since they are continually accreting from the surround- 
ing filaments (even if not through major or significant merg- 
ers). We exclude that the concentration discrepancy is caused 
by some numerical characteristics or by the physics adopted in 
the simulations. In fact, radiative simulations with no active 
and powerful feedback (such as ours) are predicted to produce 
brighter and denser cores (Borgani & Kravtsov 2011) with re- 
spect to real clusters. Therefore if any contrast in concentration 
parameter is present due to the overcooling problem affecting 
simulations, this is expected to be on the other direction. With 
the exclusion of this small difference, all other parameters agree 
with previous theoretical and observational measurements. On 
that account, we can conclude that after the extraction of the 
central part current complex theoretical models are not only 
able to fairly well reproduce standard scaling relations (Fabjan 
etaLl[20T0l |20TTj |Puchwein et 1^20081 |Battaglia et al.|2012| 



e.g.), the main properties profiles such as density, temperature, 
pressure, and X-ray surface brightness (e.g.|Nagai et al.|2007 



Fabj an et al.|2010) |Pratt et al.|2007| |Mroczkowski et al.|2 009 



Eck ert et al.|2013| Arnaud et al.|2010" 1 but also the overall ap- 
pearance of real clusters. 

This strengthen our analysis whose results can be directly 
applied to up-coming surveys. Concluding, we found that: 

• the fluctuation parameter performs better if computed 
considering the entire region within R 5 oq and a broader- 
smoothed (2') image. Also in this occasion, however, 
three disturbed halos (23% of the disturbed sample) 
fall within the regular-system range. The efficiency of 
this parameter is improved when associated with the 
asymmetry one, assuming F = 1 and A = 1.15 as divid- 
ing points within the two morphological classes; 

• the asymmetry parameter is one of the clearest indica- 
tor of morphology. The distributions of the regular and 
disturbed objects are well separated in case of A rot cal- 
culated within the R^oo region. There is no improve- 
ment in adopting variation of the asymmetry parameter 
(obtained, for example, by flipping the image instead of 
rotating it) nor in taking the maximum among differ- 
ent asymmetry parameters. This parameter is specially 
effective when it is combined with the X-ray concentra- 
tion or the third order power ratio; 

• the X-ray surface brightness concentration parameter com- 
petently divides the two morphological classes especially 
when computed using radii expressed in units of /?soq. 
The contamination is minimal, the two distributions are 
well separated as are their peaks; 

• the hardness ratio is the only parameter among those in- 
vestigated that is not able to discern among regular or 
disturbed objects. This might depend on the fact that 
both pre-merger and post-merger clusters belong to the 
last category. The thermalization states of this class, 
therefore, can be quite diverse; 

• the centroid-shift estimated in our analysis deludes our 
expectations in terms of characterizing the disturbed clus- 
ters. We believe this is due to characteristics of our sam- 
ple that present two highly complex clusters with an ex- 
tended X-ray emission symmetric to the center (CL17 
in both directions) and another object with a large sub- 
structure but still in the outskirts (CL15 both directions). 
All regular systems have a centroid-shift below 0.04 (- 
1.4 in decimal logarithm), however, also several con- 
taminants are present, including CL 15 and CL17. Im- 
provements in purity and completeness of the regular 
sample can be registered if this parameter is associated 
to others. However, even with this stratagem the com- 
pleteness of the disturbed sample remains poor (below 
50%). For the large contamination present at middle 
values of this parameter, we consider the best approach 
to identify two separate limits. For example, in our sam- 
ple, all objects with w > 0.04 are disturbed while the 
majority of systems with w < 0.01 is regular (with the 
exception of one disturbed system having w = 0.008); 

• the third order power ratio was estimated using vari- 
ous apertures. The most remarkable result was achieved 
when the parameter was computed within R^) or a fixed 
physical scale of 1000 kpc. In both cases the histograms 
related to the two distributions are completely separated. 
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The calculus of the power ratio limited to the very center 
of the cluster is inadequate for our purpose. The adop- 
tion of the peak of the power ratio profile does not im- 
prove the result substantially even if it has the gain of 
narrowing the distributions and increasing their peaks. 
The advantages of this last parameter can be enhanced 
if in the sample there are various systems with major 
substructures close to the center of the cluster. This, 
however, was not the case for the objects considered in 
this work. 

• In general terms, all morphological parameters show a 
strong correlation. For this reason, combining two or 
more parameters is a more efficient way to differentiate 
the morphological classes. 

• We defined a new morphological parameter called M. 
This is the most powerful parameter since it combines 
the strengths of all the individual parameters (it incor- 
porates detection of active dynamical history, captured 
by w, complexity and presence of substructures, typi- 
cally measured with A, F and P3). With respect to the 
individual parameter, M has two major advantages: it 
clearly has a bimodal distribution with far apart peaks 
and it is successful in distinguishing between the two 
morphological classes also when we use the median as 
reference in Eq 1, i.e. without any apriori condition on 
the boundary to be used. Studying various combination 
of the parameter M we can further determine the per- 
formances of the parameters involved. Indeed, by sim- 
ply comparing the distributions of M with and without 
the centroid-shift parameter we are able to state that this 
indicator is not suitable to uniquely categorize the dis- 
turbed class in our sample. 

Finally, it is always the case that measurements conducted 
within large scales (either or 1000 kpc) perform more ef- 
ficiently. This appear the best situation from a theoretical point 
of view, however, real clusters suffer from problematics that we 
have not considered in this paper. We refer to WeiBmann et al. 
(2013 ) for suggestions on the best approach to take in order to 
reduce the effect of Poisson noise and X-ray background. 



Appendix 

In Fig.[8]and Fig. 9] we show the samples of regular and dis- 
turbed system^] The green circles in the images are centered 
in the X-ray centroid used for rotating and flipping the images 
to compute the asymmetry parameters. Their radii correspond 
to /?5oo (see Table. 1 of R12). The cluster ID and the direction 
of the line of sight reported on the top-left corner follow the 
same nomenclature of R12. The mix class (Fig. 10 1 has been 
divided into semi-regular (including CL2 X-Y-Z, CL 4 X-Z, 
CL5 Z, CL6 X-Y-Z, CL8 X-Y-Z, CLIO X, CL1 1 Z, CL14 X- 
Y, CL16 Y, CL20 X) and semi-disturbed (including CL1 X-Y, 
CL 15 Z, CL17 Y, CL18 X-Y-Z, CL19 X-Y-Z) objects. No- 
tice that for almost all the objects, it is the case that the three 
images obtained by varying the line of sight projections belong 
to the same extended morphological class. The only exceptions 
are two clusters: CL1 and CLIO. In Table [5] we list the main 
morphological parameters for the entire set of clusters. 

7 Large resolution images are available following the link to the synthetic clus- 
ter images at http://www-personal.umich.edu/~rasia 
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Table 5 

List of parameter per each cluster evaluated with respect to R 5 qo 



CL 



M 



W 



Pi 



CL 



M 



Pi 



CL1 X 
CL1 Y 
CL1 Z 
CL2X 
CL2 Y 
CL2Z 
CL3X 
CL3 Y 
CL3Z 
CL4X 
CL4 Y 
CL4Z 
CL5X 
CL5 Y 
CL5Z 
CL6X 
CL6 Y 
CL6Z 
CL7X 
CL7 Y 
CL7Z 
CL8X 
CL8 Y 
CL8Z 
CL9X 
CL9 Y 
CL9Z 
CL10X 
CLIO Y 
CL10Z 



0.96 
0.94 
0.94 
0.73 
0.73 
0.73 
0.99 
1.00 
1.00 
1.02 
1.01 
1.00 
0.83 
0.80 
0.80 
0.80 
0.80 
0.77 
0.79 
0.83 
0.82 
0.84 
0.83 
0.83 
0.92 
0.91 
0.93 
0.80 
0.84 
0.85 



1.15 
1.11 
1.10 
0.86 
0.81 
0.86 
1.08 
1.09 
1.08 
1.14 
1.09 
1.10 
0.89 
0.82 
0.85 
0.87 
0.91 
0.88 
0.81 
0.93 
0.86 
0.92 
0.84 
0.94 
1.00 
1.08 
1.11 
0.91 
1.04 
1.32 



3.30e-07 
1.78e-07 
8.93e-08 
2.88e-07 
2.28e-07 
1.49e-07 
4.92e-08 
1.72e-07 
3.07e-08 
6.90e-07 
1.16e-07 
7.26e-07 
3.48e-08 
1.35e-07 
1.58e-07 
3.73e-07 
1.53e-07 
3.95e-08 
8.48e-08 
1.03e-07 
4.60e-08 
1.43e-07 
2.04e-07 
8.77e-09 
1.05e-09 
1.23e-07 
l.lle-07 
1.56e-07 
4.72e-06 
5.00e-06 



0.27 
0.27 
0.30 
0.41 
0.47 
0.38 
0.28 
0.28 
0.32 
0.30 
0.34 
0.31 
0.57 
0.44 
0.42 
0.39 
0.39 
0.37 
0.48 
0.42 
0.43 
0.44 
0.45 
0.48 
0.36 
0.37 
0.39 
0.39 
0.16 
0.22 



0.024 
0.024 
0.006 
0.009 
0.006 
0.008 
0.034 
0.038 
0.010 
0.034 
0.017 
0.031 
0.004 
0.021 
0.023 
0.012 
0.023 
0.019 
0.006 
0.002 
0.007 
0.008 
0.009 
0.007 
0.027 
0.008 
0.027 
0.008 
0.120 
0.104 



CL11 X 
CL11 Y 
CL11 Z 
CL12X 
CL12 Y 
CL12Z 
CL13X 
CLOY 
CL13Z 
CL14X 
CL14 Y 
CL14Z 
CL15X 
CL15 Y 
CL15Z 
CL16X 
CL16 Y 
CL16Z 
CL17X 
CL17 Y 
CL17Z 
CL18X 
CL18 Y 
CL18Z 
CL19X 
CL19 Y 
CL19Z 
CL20X 
CL20 Y 
CL20Z 



0.77 
0.79 
0.79 
1.01 
1.03 
1.04 
1.06 
1.01 
0.99 
0.72 
0.73 
0.73 
1.01 
1.01 
0.92 
0.85 
0.87 
0.86 
1.11 
1.01 
1.11 
1.27 
1.27 
1.10 
1.11 
1.06 
1.10 
0.87 
0.89 
0.88 



0.80 
0.88 
0.84 
1.34 
1.29 
1.29 
1.52 
1.44 
1.51 
0.69 
0.85 
0.79 
1.16 
1.17 
1.13 
0.82 
0.94 
0.92 
1.28 
1.14 
1.30 
1.31 
1.30 
1.34 
1.28 
1.16 
1.29 
0.94 
0.95 
0.92 



2.75e-08 
3.88e-09 
5.76e-08 
5.48e-06 
6.54e-07 
1.13e-06 
1.82e-05 
1.89e-05 
2.75e-05 
1.81e-07 
5.82e-09 
1.81e-09 
3.04e-06 
2.35e-06 
1.33e-06 
2.06e-07 
2.53e-07 
4.73e-08 
1.88e-06 
4.46e-07 
5.64e-07 
4.28e-07 
2.51e-06 
1.06e-05 
6.10e-07 
6.54e-07 
1.56e-06 
7.16e-08 
7.80e-08 
1.44e-09 



0.50 
0.54 
0.56 
0.21 
0.23 
0.23 
0.20 
0.21 
0.19 
0.48 
0.53 
0.52 
0.31 
0.30 
0.30 
0.49 
0.46 
0.45 
0.19 
0.26 
0.16 
0.31 
0.30 
0.24 
0.25 
0.33 
0.25 
0.42 
0.42 
0.43 



0.002 
0.004 
0.003 
0.064 
0.039 
0.068 
0.075 
0.080 
0.091 
0.023 
0.012 
0.016 
0.017 
0.019 
0.042 
0.003 
0.008 
0.003 
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FIG. 7. — Left panel: different combinations of the following parameters: decimal logarithm of the third order power ratio (IP3), asymmetry parameter obtained 
by rotating the image (A), X-ray surface brightness concentration (c), and decimal logarithm of the centroid shift (Iw). The corresponding (x-y) pair is written in 
each panel. Red, blue, green and magenta points denote regular, disturbed, semi-regular, and s emi -disturbed clusters. The dashed lines are the best-linear fits to 
the points (the coefficients listed in Table[5J. Right panel: distributions of the M parameter (Eq. |ll( . Red and blue histograms refer only to regular and disturbed 
systems while red-magenta and blue-green distributions include semi-regular and semi-disturbed objects, respectively. 
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FIG. 9. — X-ray soft images (0.7-2) keV of the cluster visually classified as disturbed. The green circle indicating S500 has been center on the X-ray centroid. 
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